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CONVEX HULLS OF ORBITS AND ORIENTATIONS OF A MOVING 

PROTEIN DOMAIN 
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' Abstract. Wc study the facial structure and Caratheodory number of the convex hull 

of an orbit of the group of rotations in acting on the space of pairs of anisotropic 
symmetric 3x3 tensors. This is motivated by the problem of determining the structure 
of some proteins in aqueous solution. 

m 

o 

Introduction 

^ . The most aesthetically appealing polytopes arise as convex hulls of orbits of finite groups 

, acting on a vector space. These include the platonic and archimedean solids and their 

higher- dimensional generalizations, such as the regular polytopes [7]. In contrast, the 
analogous objects for compact Lie groups have not attracted much study. We investigate 
^ . convex hulls of orbits of the group 5*0(3) in a particular 10-dimensional representation. 

^ [ Our motivation comes from an algorithm to understand the fold of some proteins. 

Certain proteins, such as calmodulin [3], consist of two rigid domains connected via 



■ a region that is flexible in aqueous solution (i.e. under physiological conditions), and 

^ , the problem is to determine the relative position and orientation of these two domains. 

^ I Calmodulin, as many other proteins, incorporates metal ions into its structure. When 

' a paramagnetic ion is substituted, it interacts with the magnetic field of dipoles within 

O ■ the protein via its magnetic susceptibility tensor x- Part of this interaction, the residual 

dipolar coupling, may be inferred from nuclear magnetic resonance data and depends 
' solely upon the relative orientation of the two domains. When the relative orientation of 

^ ■ the two domains is not constant, we infer the mean magnetic susceptibility tensor x from 

these data. 

We model this relative orientation by a probability measure p on the group 5*0(3) of 
rotations of M^. Then x is the average with respect to p of rotations of x- Recovering p 
from X is an ill-posed inverse problem. Nevertheless, x contains useful information about 
p. Gardner, Longinetti, and Sgheri [TU] gave an algorithm to determine the maximum 
probability of a given relative orientation of the two domains. Since x lies in the convex 
hull V of the orbit of x under the group of rotations of M^, it admits a representation 
X = T^jVjRj-X^ where the sum is finite, J2jPj — ^ with pj > 0, Rj is a rotation in M^, 
and Rj-X is the action of Rj on the tensor x- The minimal number of summands needed 
to represent any point x in F is the Caratheodory number of V. 
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It is often possible to substitute a different metal ion into the protein with a different 
susceptibility tensor x'. Repeating the measurements gives a second mean tensor x' which 
is the average of rotations of x' with respect to the measure p. Combining this with x 
gives more information about p. Longinetti, Luchinat, Parigi, and Sgheri [12] adapt the 
algorithm of [10\ when there are two or more metal ions and show how this can be used 
to better understand the structure of calmodulin. Their algorithm uses some knowledge 
of the convex hull V^'"^ of the orbit of the pair {x, x') under the group of rotations. 

We study the Caratheodory number and facets of V^'"^. When x and x' linearly 
independent, V has dimension 10. We call a subgroup of SO (3) which stabilizes a line in 
a coaxial group, and a face of V^'"^ which is stabilized by a such a subgroup a coaxial 
face. Our main result is the following. 

Theorem 15.61 Faces ofV^'"^ have dimension at most 6. The coaxial faces of V^'^ form 
a 3 -dimensional family whose union is a 9-dimensional subset of the boundary of V^'"^ if 
and only if x and x' have distinct eigenvectors. In that case, almost all coaxial faces have 
dimension 6, have Caratheodory number 4, and are facets of V^'"^. 

Our main result implies that the Caratheodory number of V^'"^ is at most 8. This is an 
advance over [10], where it has been bounded between 4 and 10 inclusive. 

We are unable to show that the boundary of V^'"^ is the union of coaxial faces when x 
and x' have distinct eigenvectors, but conjecture that this is the case. As a consequence 
of our main result, we also conjecture that the Caratheodory number of V^'^ is at most 5. 

Given any number of tensors (xi, • • • ,Xn) we may define the convex hull V^''"'^ . 
In Subsection 13.21 we prove that dimV^'"'^ is 5 times the dimension of the span of 
(Xi, • • • , Xa^)- Iu the text, we will omit the superscripts from our notation for the convex 
hull. 

Magnetic susceptibility tensors are 3x3 symmetric trace zero matrices and form a 
5-dimensional irreducible representation of the group S0{3) of rotations in M^. More 
generally, one could study the convex hulls of orbits of compact groups. We were surprised 
to find that very little is known about such convex bodies, particularly their Caratheodory 
numbers and facets. We hope that our work will stimulate a more thorough study of 
convex hulls of orbits of compact groups. 

In Section [T], we describe the motivation for this work from protein structure. In 
Section [2] we discuss group actions and in Section [3] convex hulls of orbits. In Section H] 
we complete the analysis of [TU] in the case of one metal ion. In Section [S] we analyze the 
case of two metal ions and deduce Theorem 15.61 

1. Application to protein structure 

Proteins are large biological molecules synthesized by living organisms. The genome 
of an organism contains the chemical formulae for its proteins. Currently, hundreds of 
organisms (including man) have had their genomes mapped, and such chemical formulae 
are readily available. An important step towards inferring the biological function of a 
protein from its chemical formula is to determine its 3-dimensional structure, or its fold. 

About one third of all proteins incorporate metal ions into their structures. The fold of 
these proteins may be inferred from nuclear magnetic resonance, which can measure the 
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interactions between paramagnetic metal ions and dipoles within the protein. The main 
quantities that can be measured are the pseudo contact shifts (PCS) pQ and the residual 
dipolar coupling (RDC) [17]. In this paper we only deal with the RDC 

The residual dipolar coupling between a paramagnetic ion and a dipole formed by atoms 
a and b within the protein depends upon the vector displacement r from the atom a to 
the atom b and the magnetic susceptibility tensor x of the metal ion, which is a 3 x 3 
symmetric matrix. The RDC interaction has the following vector formula 



:i.ii 



c 



T 

-r xr 



C 



3||r||3 



Trace (x) . 



Here, C is a constant and ||r|| is the length of the vector r. This only depends upon the 
relative orientation of the dipole and metal ion, and so the RDC data may be used to 
infer this relative orientation. 

Writing X = Xo + |Trace(x)/3, where is the 3x3 identity matrix and Xo is the 
trace-free or anisotropic part of Xi this formula becomes 

X C rp 

We assume henceforth that x = Xo is anisotropic. 

The fold of the protein is usually unique, in the sense that small variations of the shape 
are allowed. There are proteins which however exhibit large variations of shape under 
particular conditions. A widely-studied example is calmodulin, which has two rigid do- 
mains, called the N-terminal and C-terminal domains, connected by a short flexible linker. 
The N- and C-terminal domains are assumed to be rigid bodies with known structures. 
Figure [H obtained with Molmol [11] , shows calmodulin in two different orientations. 

N-terminal 



C-terminal 




linker zone 



Figure 1. Two orientations of calmodulin 



The binding site of the metal ion in calmodulin belongs to the N-terminal domain. 
The measured RDC of pairs of atoms belonging to the N-terminal domain can be used 
to obtain a good estimate of x- The measured RDC of pairs of atoms belonging to the 
C-terminal domain can be used to study the relative orientation of the two domains. 

Let us model the relative orientation of the N- and C-terminal domains with a rotation 
R. Then there is an unknown probability measure p on the set 5*0(3) of rotations such 
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that the mean RDC 6 of the pair of atoms a, h in the C-terminal domain is given by 

(1.2) J = ^ [ {Rrfx{Rr)dp{R) = /^^r^X^ , 
where the mean magnetic susceptibility tensor x is 

(1.3) X = f R^xRdp{R). 

JSOz 

This tensor x can be estimated from the RDC of several dipole pairs in the C-terminal 
domain. The experimental measures show that in terms of difference of eigenvalues, x is 
between 5 and 20 times smaller than x 0- This indicates that p is not a point mass, that 
is, the C-terminal domain moves with respect to the N-terminal domain. 

The availability of N distinct mean susceptibility tensors Xk with respect to different 
metal ions /c = 1, . . . , iV increases the information about p for up to 5, see for instance 
|14j and [121 Theorem 3.2]. However, even the exact knowledge of 5 mean tensors Xk (i-^- 
25 real numbers) does not allow the exact reconstruction of the probability measure p. 

An approach to extract information from the mean tensors is to define Pmax(-R) as the 
maximal fraction of time that the C-terminal can stay in a particular orientation i?, yet 
still produce the measured mean tensors. Orientations with a large Pmax agree with with 
what are thought to be the most favored orientations of the C-terminal |12j . 

The calculation of Pmax can be performed geometrically if only RDC is considered [12] . 
In the combined PCS+RDC case more information is added [1], however the calculation 
of Pmax can only be performed numerically. For the convergence and efficiency of the 
algorithm, the minimal number of orientations needed to reconstruct any admissible set 
of mean tensors Xk should be used. Experience suggests that adding the PCS data does 
not increase the actual number of orientations needed [1], so the Caratheodory number 
for the RDC case may be used as a basis for the numerical minimization. 

2. Groups actions and anisotropic tensors 

We first recall some basics about representations of compact groups, then consider the 
action of the group 5*0(3) of rotations in M? on the 5-dimensional space of anisotropic 
3x3 symmetric tensors, and finally investigate the coaxial subgroups of 5*0(3). 

2.1. Representations of compact groups. This material may be found in the book 
of Brocker and tom Dieck [5l Ch. II]. Let G be a compact group, such as 5*0(3). A 
representation of G is a finite-dimensional vector space W on which G acts by linear 
transformations. That is, we have a group homomorphism p: G GL{W), where GL{W) 
is the group of invertible linear transformations on W. For g E G and w G W , write g.w 
for p{g){w). 

A representation 14^ of G is irreducible if its only G-invariant subspaces are {0} and 
W . Every representation of G decomposes as the direct sum of irreducible representations 
which is unique in the following way. Given a representation W oiG and a positive integer 
Z, let be the /-fold direct sum of W , 

^ ^ ^ 
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Elements g oiG act diagonally on elements w = {wi, . . . , wi) of W^', g.w = {g.wi, . . . , g.wi). 
Suppose that Wi, W2, ... is the list of irreducible representations of G. If ?7 is a represen- 
tation of G then there exist unique integers /i, /2, . . . such that 

as representations of G. If f/j is the subrepresentation of U mapped to the summand 
under this isomorphism, then Ui does not depend on any choices and is called the 
isotypical component of U corresponding to Wi. If /j > 0, then we say that U contains 
Wi- Furthermore, if f/' C f/ is a subrepresentation, then the ith isotypical component of 
U' is U' n Ui, which is also the image of U' under the projection to Ui. 

Haar measure is a G-invariant measure fi on G with 1 = J^d^lg). Given a linear 
function L : W ^ M., where is a representation of G, we may average L over orbits of 
G to get a new linear function L', defined by 

L'{x) := / L{g.x)dfi{g) . 
Jg 

Since L' is constant on orbits of G, it is G-invariant. This association L L' is called 
the Reynolds operator. It is an important tool for analyzing G-representations. 

Another key tool is Schur's lemma. A linear map (f.W^U between representations 
of G is a G-map if for aX\ w & W and g & G, we have g.(p{w) = (p{g.w). Let HomG(Vr, U) 
be the space of G-maps. A division algebra is a finite-dimensional associative algebra in 
which every non-zero element is invertible. 

Schur's Lemma. If W U are irreducible representations of G, then HomdW, U) = 
and HomciW, W) is a division algebra which contains M. 

Proof. Let (f: W ^ U he a G-map. Then the kernel of </9 is a subrepresentation of W and 
so it is either or W, and the image of if, which is a subrepresentation of U, is either 
or U. Examining the possibilities leads to the conclusions. □ 

There are exactly 3 division algebras which contain M: The real numbers M, the complex 
numbers C, and the quaternions H. An irreducible representation of G has real, 
complex, or quaternionic type, depending on EndG(Vr) := B.omG(W,W). 

Example 2.1. Consider the group 5*0(2) of rotations of M^, 

Ends'o(2)(K^) consists of those 2x2 matrices M such that MRg = RqM, and so 

End50(2)(M^) = m(J J) + m("J J ) . 

This is isomorphic to C (we send (^"^ to V— 1), so this representation of S0{2) has 

complex type. It is the defining representation Ui of 5*0(2). If we identify M? with C and 
5*0(2) with the circle group := {e*^ | < ^ < 27r}, then the action is scalar multipli- 
cation by elements of 5*^. For any positive integer /c > 0, let Uk be the representation of 
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on C (identified with M?) where z E acts as multiphcation by z^. These all have 
complex type. 

Example 2.2. The orthogonal group 0(2) contains 50(2) as well as the cosets of reflec- 
tions 

S0{2).(^1 "J) = {Ro-(^~1 J) I ^e[0,27r)| . 
The defining representation Ui of 0(2) on has real type, as 

(oi)'(i o)"(io)^(-i o)"(i o)'(oi)' 

and so Endo(2)(K^) = ^ " (o J) - ^• 

In the trivial representation f/g = M of 0(2), elements act as multiplication by 1. For 
a positive integer k, define the map (pk'. 0(2) 0(2) C GL(2,]R) by (pk^Re) = Rue and 
'^^("o i) ~ (~o i)' '^^^^ defines the representation Uk of 0(2), which has real type. 
Restricting to 5*0(2) gives its representation Uk of complex type. 

2.2. Rotations of anisotropic tensors. Let ei = (1,0,0)"^, 62 = (0,1,0)"^, and 63 = 
(0, 0, 1)^ be the standard basis of M'^. 

The special orthogonal group 50(3) is the group of rotations in R^. It consists of 3 x 3 
real orthogonal matrices with determinant 1, 

50(3) := [ReW^""^ \RR^ = 1 and deti? = l}. 

Let R G 50(3), and let Tr50(3) be the tangent space to 50(3) at the matrix R. Let I 
be the identity matrix, then T/50(3) is the space SO3 of skew symmetric 3x3 matrices, 
which is the Lie algebra of 50(3). That is, 

Tj50(3) = / + SO3. 

Elements R G 50(3) act on 3 x 3 symmetric matrices (tensors) x by conjugation, R.x '■ = 
RxR^^- This preserves the trace of and so 50(3) acts on the space W of anisotropic 
(trace-zero) tensors, a 5-dimensional irreducible real representation. We introduce some 
useful coordinates for W. A point (f , w, x, y, z) G corresponds to 

( V 0\ /OwxX /OO 0\ 
(2.3) x{V:W,x,y,z) := -| + U; + y z \ . 

\0 0-|/ \a;00/ \Q z -y j 

Observe that ei is an eigenvector for x{v,w,x,y, z) if and only if u; = x = 0, 62 is an 
eigenvector if and only if w = z = 0, and 63 is an eigenvector if and only if x = z = 0. 



tXhis left action (if i?, 5* G S0{3), then R{S.x) — RS-x) is equivalent to the implied action in (|1.3p . 
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2.3. Coaxial subgroups. The coaxial subgroup is the set of rotations fixing a line 
in with direction e. It is isomorphic to the orthogonal group 0(2). Its identity 
component Q'^ is isomorphic to SO (2) and consists of rotations about the axis e, while 
the other component Q~ consists of refiections in axes orthogonal to e. For example, let 

/lO 0\ /-100\ 

(2.4) i?ei,0 := cos^ -sin^ and Re,,^ := 0-10. 

\0 sin^ cos^y ' \ 1 / 

If we fix / perpendicular to e and let 9 run over all angles, then i?e,e and RefiRf,-K give all 
elements of Qe- 

We consider the action of a coaxial subgroup Qe on W . For this, suppose that e = ei 
and let Re^fi act on ^(w, w, x, y, z). This gives the tensor xW ^ '^'y ^'y v'^ where 

v' = V , 

hi = f = Re and 

\x J y smt) COS& J \X J V-^/ 

fy'\ _ /cos 2^ -sin2^\/y\ _ f y\ 

\z'J ~ \sm2e cos 26 J {z J ~ ^^'{zj- 

Thus i?ei,6» acts trivially on the coordinate v, by rotation through the angle 6 on the 
vector {w, x)'^, and by rotation through 26 on the vector {y, z)'^ . Note that Re^^-n sends 
x{v,w,x,y, z) to x(^) w,—x, y,—z). Thus, if we restrict the action of 5*0(3) on W to 
Qei — 0(2), then it decomposes as a sum of irreducible representations 

(2.5) W = Uo ® Ui ® U2. 

This decomposition corresponds to the coordinates ( 12. 3p . Projection to the trivial submod- 
ule f/o = ffi is, up to a scalar multiple, the unique Qe-invariant linear function L : — > M. 

3. Convex hulls of orbits 

Let O be an orbit of a compact group G in a representation W of G. The convex hull 
V of O is all points of W which are convex combinations of elements of O, 

XlVi + X2V2 H h XnVn , 

where G O, and the non- negative numbers Aj have sum 1. The set ^ is a 

compact convex set, hence a convex body. 

3.1. Faces and Caratheodory number of V. A face F of ^ is the subset of its bound- 
ary where some linear function L achieves its maximum on V, 

F := {v eV \ L{v) > L{u) for all ueV} . 

We say that L supports F and also that the hyperplane L{x) = L{F) supports F. (Here, 
L{F) is a constant.) The tangent spaces to O of its points lying in F are contained in 
any hyperplane supporting F. Such a tangent space TyO at the point v is 

V + Q.V, 

where the action of the Lie algebra g is the derivative of the action of G. The face F is 
proper if F V. When V is full-dimensional so that dim V = dim W, this is equivalent 
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to L 7^ 0. A facet is a maximal face and a vertex is a minimal face. Vertices are not a 
convex combination of other points of V. 

Lemma 3.1. The vertices of V are exactly the points of O . 

Proof. The vertices of V are a subset of O := G.x. Let g.x G O and suppose that it is a 
convex combination of vertices, 

g.x = Xigi.x + X2g2-x ^ hA„fif„.x. 

Multiplying by gig~^ expresses the vertex gi.x as a convex combination of points of O. 
Thus n = 1 and g.x = gi.x is a vertex. □ 

When V has dimension d, Caratheodory's Theorem [6], see e.g. [T6l Theorem 1.1.4], 
implies that any point x of is a convex combination of at most d+1 vertices. The 
Caratheodory number of V is the minimum number n such that any point x G is a convex 
combination of at most n vertices. For example, a ball in has Caratheodory number 2, 
while a (i-simplex has Caratheodory number d + 1. Fenchel [9], see e.g. [151 Theorem 1.4], 
showed that the Caratheodory number is at most d when the set of vertices is connected. 
More useful for us is a recursive bound, which is immediate from the observation that any 
point of V is the convex combination of any vertex and some boundary point. 

Lemma 3.2. The Caratheodory number of a convex body V is at most one more than the 
maximal Caratheodory number of its facets. 

Suppose that S' is a closed (hence compact) subgroup of G which stabilizes a face F of 
V, that is s.F = F for all s e S. 

Lemma 3.3. When F is proper, there is a non-zero S -invariant linear function on W , 
and W contains the trivial representation of S. 

Proof. Let L : — M be any linear function supporting F with L{F) = £. Let L' be the 
image of L under the Reynolds operator for S. For m G -F, we have L'{u) = £, as F is 
S'-stable and L'{u) is the average of L over the orbit of S through u. 

Suppose that w ^ V \ F. Then L{w) < L{u) = £ and L{s.w) is bounded away from 
£ as S.w is compact and disjoint from F. In particular, this implies that L'{w) < L'{u), 
which shows that the S'-invariant linear function L' supports F and that L' ^ 0. Such an 
5'-invariant linear function must factor through the trivial isotypical component of W as 
a representation of S. This completes the proof. □ 

3.2. The dimension of V. If w = {wi, . . . , Wi) G W\ then we write d{w) for the dimen- 
sion of the linear span of the components Wi, ... ,wi of w in W. 

Lemma 3.4. Suppose that W is an irreducible representation of a group G having real 
type. If w E W\ then the linear span U of the orbit G.w in is isomorphic to W'''^'^'^ . 

Proof. Write k := d{w). We may assume that wi, . . . ,Wk form a basis for the linear span 
of wi, . . . ,wi. Let A = (aij) G Mat;xfc(K) be the matrix which writes the components 
of w in terms of this basis, Wi = for z = 1, . . . , /. For each i = 1, . . . , Z, let 
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(fi: U ^ W he the projection to the ith coordinate. Since ipi{g.w) = g.Wi, we have 

k 

(3.5) yj. = ^ aijLfj for i = l,...,l. 

i=i 

This matrix A defines a G-map A: by 

(3.6) A ■ {wr,...,Wk) ^ [Y, 

\ j j j 

Composing the map i/j := {(pi, . . . , ipk) '■ U — ^ W'' with A: gives the identity 

map on f/: By (13.51) . for w E U, we have 

We show that the map (13. 6p is injective and thus is an isomorphism. A hnear map 
L: — > M induces maps L'^: W'' and : M', which commute with A. If 

7^ w G W^, then there is some hnear map L: W with L^{w) ^ 0. Since A has full 
rank k, A{L''{w)) ^ 0. But this implies that A{w) ^ 0, as A{L''{w)) = L\A{w)). □ 

Lemma 3.7. Suppose that W = Wl'^ © ■ ■ ■ © is the decomposition of a representation 
W of G into isotypical pieces, each of which has real type. Let w = {wi, . . . ,Wm) G W , 
where Wi is the component of w in W^^ . Then the dimension of the convex hull V of the 
orbit G.w is 

d{wi) ■ dimiyj . 

IfW does not contain the trivial representation, then lies in V. 

Proof. If W contains the trivial representation, assume that it is Wi. Since g.{w — Wi) = 
g.w — Wi, we see that the orbits G.w and G.{w — Wi) are isomorphic, and the same it true 
for their convex hulls. Thus it it no loss to suppose that wi = 0, which is equivalent to 
assuming that W does not contain the trivial representation. 

The linear span U of the orbit G.w is the direct sum of its projections to the isotypical 
components W^' of W. Each projection is the linear span of G.Wi, which by Lemma [3.41 
is isomorphic to wf^^^\ Thus U has dimension ^id{wi) dimPVj. We may replace W by 
this linear span, and therefore assume that the orbit G.w spans W. 

The convex hull of G.w lacks full dimension only if it lies in some hyperplane H not 
containing the origin. Suppose that this is not the case and let B be the convex hull of 
G.w and the origin. Then is a proper G-stable face of B and so by Lemma 13.31 W 
contains the trivial representation, which is a contradiction. 

If ^ V^, then there is some linear function L which is bounded above on V. But 
then the image L' of L under the Reynolds operator is non-zero on V. This implies that 
L' 7^ 0, and so W contains the trivial representation, a contradiction. □ 

Example 3.8. Lemmas 13.41 and 13.71 do not hold if the representation W has complex 
type. For example, let G = S0{2) and W = Ul with k,l > 1. Identifying W with 
and SO (2) with the circle group S^, elements z E S^ act on as scalar multiplication 
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by z'^. Thus the hnear span S0{2).w for w & Ul is a. complex hne, and therefore has real 
dimension 2, and not 4 as Lemma [3.41 predicts for general w & W when l,k > 1. 

In particular, if = f/i © f/|, and w & W is general, then the linear span of G.w has 
complex dimension 2 and thus real dimension 4. 

4. One metal ion 

Let W be the space of symmetric 3x3 anisotropic tensors, a 5-dimensional irreducible 
representation of 5*0(3) of real type. For each unit vector e G M^, there is a linear function 

(4.1) Le : W3X ^ {e,xe) = eV^M. 

If e is an eigenvector of then Lf,{x) is its eigenvalue. In general, i^e(x) li^s between 
the maximum and minimum eigenvalues of x- Note that Lg is Qe-invariant. By the 
decomposition (12.51) of W into irreducible Qe representations, any Qe-invariant linear 
function is a scalar multiple of Lg. 

As a matrix, a tensor in W has real eigenvalues and its eigenvectors form an orthonormal 
basis for M.^. Fix a non-zero anisotropic tensor x ^ with maximum eigenvalue M > 
and minimum eigenvalue m < 0. The intermediate eigenvalue of x is — M — m, and we 
have —^>m> —2M. The orbit O of x under SO (3) consists of the anisotropic tensors 
with maximal eigenvalue M and minimal eigenvalue m. It is a manifold whose dimension 
we determine. 

Proposition 4.2. The orbit is ?>- dimensional unless x has an eigenvalue of multiplicity 
2 and then it is two-dimensional. 

Proof. The dimension of O is equal to the dimension of any of its tangent spaces. Since 

O = [RxR^ I R e 50(3)} , 
the tangent space T^C at a point ^ G C is the affine space 

(4.3) i + {ri + ir^\re SO3} . 

Indeed, consider the action in ^ of an element J+r of the tangent space T/50(3) = /+SO3: 

(/ + r)^(J + rY = ^ + r^ + ^r^ + r^r'^ . 

Discarding the term which is quadratic in SO3 gives ( 14. 3p . 

It suffices to determine the tangent space to O at the point x- We may suppose that x 
is diagonal, and let r be a general element of 503, 

/M 0\ /O -a -h\ 

X = I — M — m and r = I a — c , 

\0 mj \b c J 

where a,b,c e R. Let a := M + p := M - m, and 7 := -m - f^. Then /3 > and 
a, 7 > with a = only when the eigenvalue m has multiplicity 2 and 7 = only when 
the eigenvalue M has multiplicity 2. We see that T^O is the affine subspace of W 

/ 2aa 6/3 \ 

(4.4) ^ + \2aa 2c7 , 

\bf3 2c7 / 
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where a,b,c G M. This is 3- dimensional unless either a = or 7 = 0. □ 

Let V be the convex hull of the orbit (9 of x G W. By Lemma [3^ this is a 5- dimensional 
convex body. 

Lemma 4.5. If x ^ ^ csnc? e G M'^ ^s a unit vector, then we have 

(4.6) m < Leix) < M. 

In fact, V is the set of symmetric anisotropic tensors satisfying fl4.6p flUl Theorem 3.3]. 
Proof. Le{x) lies between the maximum and minimum eigenvalues of x- Thus the in- 
equality (14. 6 p holds for x iii the orbit of x- Since a general element of is a convex 
combination of tensors in the orbit of we deduce (14.61) . □ 

A coaxial face of is a face which is stabilized by some coaxial subgroup Qe- By 
Lemma 13.31 a coaxial face stabilized by Qe is supported by a non-trivial Qe-invariant 
linear function. As we noted earlier, this linear function is necessarily a scalar multiple of 
Lg. By the inequality (14. 6p . there are two possibilities for such a coaxial face, 

:= {xeV\L,{x) = M} and 

(4.7) 

:= {xeV\L,ix) = m} . 

The coaxial face F^^ consists of tensors x G having e as an eigenvector with eigenvalue 
M and tensors in have e as an eigenvector with eigenvalue m. 

We now describe the facets of V and determine its Caratheodory number. As in the 
proof of Proposition 14. 2[ set a := M + ^ > and 7 := — m — ^ > 0. 

Lemma 4.8. The boundary of V is the union of coaxial faces (14. 7p where e ranges over 
all unit vectors in M.^ . A nonempty intersection of two coaxial faces lies in the orbit of 
SO (3). Each face F^^ is a circle of radius 7 and each face F^ is a circle of radius a. 
When X has a repeated eigenvalue so that either a or vanishes, then the corresponding 
coaxial face degenerates to a point. 

A consequence of Lemma 14.81 is that the coaxial faces are maximal faces, and are 
therefore facets. 

Proof. By Theorem 3.3 of [10], V is the set of anisotropic tensors whose eigenvalues lie 
in the interval [m, M], and so its boundary consists of tensors x either having maximal 
eigenvalue M or having minimal eigenvalue m. This shows that the boundary of V consists 
of coaxial faces, which are thus the facets of V. 

We show that the intersection of two coaxial faces lies in the orbit of x- Suppose 
that X lies on two different coaxial faces. If these are F^^ and FJ^, then e and / are 
eigenvectors of x with eigenvalues M and m, respectively. The third eigenvalue of x is 
—M — m, and so x lies in the orbit of x- If the two faces have the form F^'^ and Fj^ with 
e and / linearly independent, then the eigenvalue M of x has multiplicity 2 and its third 
(smallest) eigenvalue is — 2M. Since x ^ ^) this smallest eigenvalue is bounded below 
by m; as m > — 2M, we see that m = —2M and so again x lies in the orbit of x- The 
argument is similar if the two faces are F^ and F'p. 

The coaxial face F^^ of V consists of tensors x G having e as an eigenvector with 
eigenvalue M. Since each point of the boundary of F^ lies in some other coaxial face. 
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this boundary lies in the orbit of x ^i^nd is necessary an orbit of Qe- We need only 
consider the case when x ^ so that this boundary is Qe-X- Suppose that e = ei and 
X = 0, 0, 7, 0), in the coordinates (12.31) . Here, 7 = — ^ — rn>Q. As in Section [^75| 
elements of Qe act on x by rotation of the vector (7, 0) formed by the last two coordinates, 
and thus F^^ is a circle of radius 7, which degenerates to a point if 7 = 0. 

We omit the similar arguments for F^. □ 

Theorem 4.9. If zero is not an eigenvalue of x, then V has Caratheodory number?), and 
when zero is an eigenvalue, V has Caratheodory number 2. 



When zero is not an eigenvalue of this is the main result about V from jlOj . 

Proof. First suppose that zero is not an eigenvalue of x- Every facet has Caratheodory 
number 2, as it is a circle. So by Lemma [3.21 V has Caratheodory number either 2 or 3. 
By Lemma 13. 7[ lies in V . If V has Caratheodory number 2 then there exist A G [0, 1] 
and R,S e 50(3) with = \R.x + (1 - A)^.^. Multiplying by this becomes 

= Ax + (1 — A)i?.x, for a different rotation R G S0{?>), and so 

(4.10) -Ax = (1-A)i?.x. 

Suppose that x is diagonal. Then (14.101) implies that R.x is also diagonal. 

If M 7^ —m so that is not an eigenvalue of x, then one of the diagonal matrices —Ax 
and (1 — \)R.x has two positive entries and the other has two negative entries, which is 
a contradiction. Thus if is not an eigenvalue of Xi then V has Caratheodory number 3. 

Now we assume that is an eigenvalue of x- We will show that the image of the map 
[0, 1] X S0(3) X 50(3) ^ V which takes (A, R, S) to \R.x + (1 - \)S.x meets each 50(3)- 
orbit in V and is therefore surjective. They key point is that two tensors are in the same 
orbit if and only if they have the same characteristic polynomial. 

The characteristic polynomial of a trace- zero matrix x with eigenvalues s, t, — s — t is 

- x{st + t^ + s^) + {sH + st"^) . 

The constant term is — det(x), while the coefficient —a of x is the sum of the pairwise 
products of eigenvalues, which is an invariant of the matrix. 

Scaling x, we may assume that its eigenvalues are 1,0, and —1, so that V consists of 
tensors x with eigenvalues s, t, — s — t lying in the interval [—1, 1]. The set of such 
pairs (s,t) are the points of the hexagon of Figure [2J The three lines through the origin 
t = s, t = —s/2, and t = 2s divide the hexagon into six quadrilaterals and permutations 
of the eigenvalues permute these quadrilaterals. We leave it to the reader to check that 
{s,t) I— > (a;,det) is a one-to-one mapping of each quadrilateral onto the region shown in 
Figure [2], which is 

(4.11) X := {(a,det) G | 27det^ < 4a^ and a < 1 - | det | } , 

and is bounded by the curves det = 1 — a, det = a — 1, and 27det^ = 4a^. 
Consider matrices of the form x(A, 0, r) := \R{6).x + (1 — A)S'(r).x, where 



m 



cos 9 sin 9 
— sin 9 cos 9 
1 



and S{t) 



cos r sin r 

1 
- sin r cos r 
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t = ~2s 




s 
2 




Figure 2. Eigenvalues and invariants of tensors in V. 

The invariants (a, det) of 9, r) are 

(1- A(l-A)(4sinV + sin2^-2sin2rsin2e), A(l - A) sin^ ^(1 - 2A sin^ r)) . 

If we let u = sin^ 9 and v = sin^ r, then the set of invariants of xl-^? 9, r) for all (A, 9, r) 
are the image of the unit cube [0, 1]^ under the map 

/: i\u,v) I — ^ (1- \){'iv + u-2uv), A(l - A)m(1 -2At;)) . 



We show that the image of / includes that part of X (14.111) where det > 0. This will 
complete the proof, as replacing x by — x = S{^).x E V in our definition of / changes 
the sign of the determinant and does not change the invariant a. 

Figure [3] shows subsets of the faces v = 1 and m = 1 of the cube which include the 




Face V 




Face u 



Figure 3. Subsets of faces of the cube 



segments A = 1/3 and have boundaries the indicated curves. The map / is one-to-one on 
the interior of each region, and the images cover that part of X with det > 0, meeting 

only along the curve det 

line 3 det 



(a 



\){(y — ly for I < a < 1, as shown in Figure HI The 

, which is tangent to the boundary curve 27det^ = Aa^ at the point 
(|, ^), is the image of the lines A = | in Figure [31 This completes the proof. □ 

Lastly, we identify the hyperplanes supporting facets of V. By Lemma 14.81 the faces 
and are 2-dimensional, unless x has an eigenvalue with multiplicity 2, and in that 
case exactly one face is degenerate. 
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det = ^A(«-i)(«-l)2 

Figure 4. Image of subsets of Figure [3] 

Proposition 4.12. The facets and have a unique supporting linear function, 
unless they are degenerate. 

Proof. It suffices to determine the hyperplanes supporting faces which contain x- Any 
hyperplane supporting the vertex x contains the the tangent space T^C at x to the orbit 
O through X- Choose coordinates so that x = x(^) 0, 0, 7, 0) is diagonal (7 = — ^ — m) 
so that X e n F^. Recall that T^O (glD is the affine 3-plane in W 

X(M,*,*,7,*) , 

where * represents an arbitrary real number. If (xi,...,X5) are the coordinates (12 .Sp . 
then hyperplanes containing T^C have equation J2i (^i^i = c, where 

(4.13) C2 = C3 = C5 = 0, and ciM + C47 = c. 

If 7 7^ 0, then F^^ is nondegenerate and contains the additional point 

/MO \ 

i?e„f = x(M,0,0,0,7) = -f 7 , 

VO 7 -fj 

which imposes the further condition CiM = c on a support hyperplane to F^^ . Thus 
Ci = c/M and C4 = 0. Setting c = M so that Ci = 1, we see that the support hyperplane 
to Fg^^ is defined by Xi = M, which is ivei(x) = 

If u = M + Y ^ 0, then F^ is nondegenerate and it contains the point 

-f ^ o\ 

^ -f = x(-f,i^, 0,-^,0), 

mj 

and so a support hyperplane to F^ must satisfy (14.131) and also — Ciy — 04^ = d. 
Subtracting these equations and dividing by u, we see that ci = 2c4, and so 

Xi + 2^4 = —2m , 

is the support hyperplane to F^. Note that xi + 2x4 is — 2Le3(x)- D 
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5. Two Metals 

Let Xi) X2 £ W he linearly independent anisotropic tensors and set x := (xi? X2) G W^. 
By Lemma [3.71 the convex hull V of the orbit O := S0{3).x is a 2 ■ 5 = 10- dimensional 
convex body containing the origin. Its boundary is 9-dimensional. 

We study the facial structure ofVc. W^. One tool will be a family of S'0(3)-equivariant 
maps tTq,: W"^ — > W. We first determine the dimension of the orbit, show that the 
maximum dimension of a facet is 6, and then define coaxial faces. Our main result is that 
coaxial faces are facets if Xi and X2 have distinct eigenvectors. In that case, almost all 
coaxial faces have dimension 6 and Caratheodory number 4. We are unable to rule out 
the existence of other facets, but we conjecture that there are no other facets. 

Let Span(x) G W he the 2-dimensional subspace of W spanned by xi X2- The 
structure of V depends only on Span(x). Indeed, if XijX2 ^ Span(x) are linearly inde- 
pendent, then there is a 2 x 2 invertible matrix A = (ajj) such that 

xi = fliiXi + 0.12X2 X2 = 021X1 + ^22X2- 
This induces an S'0(3)-isomorphism W^: 

3 {wi,W2) I — > {aiiWi + ai2W2, 021^1 + 022^2) e W"^ 

which sends V to the convex hull of the orbit of (x'i5X2)- This is nothing more than a 
change of coordinates on W^. 

Any non-zero vector a = (ai, 02) ^ gives an S'0(3)-map 

(5.1) TT^ : — > W 

defined by 7!'a{wi,W2) := aiWi + a2W2- Write Wa for iTaiw). In particular, Xa '■= ^^aix) ^ 
Span(x). Set Va G W to he the convex hull of the orbit Oa '■= S0{3).Xa- Since 
7ia{SO{3).x) = 50(3). Xa, we have TTaiO^) = Oa and Va = HaiV). We compute the 
dimension of the orbit O^. 

Theorem 5.2. dimO^ = 3. 

Proof. We will show that dim Oa = 3 for some a G M^. As iVaiOy^) = Oa, this implies 
dim0^ > 3. Since dim 5*0 (3) = 3, we have dim (9^ < 3 and so dim (9^ = 3. 

By Proposition \A.2\ the dimension of Oa is 3 if and only if Xa has distinct eigenvalues. 
If either xi or X25 say xi; has distinct eigenvalues, then dim(9(i o) = 3 and we are done. 
Suppose the contrary, that neither Xi nor X2 has distinct eigenvalues. That is, for each 
i = 1, 2, Xi has a 2-dimensional eigenspace with eigenvalue ctj. Since cannot be a repeated 
eigenvalue, neither ai nor a2 is zero. These eigenspaces must meet, so Xi and X2 share an 
eigenvector, which is an eigenvector for the nonzero tensor Xa '■= «2Xi ~ ^1X2 ^ Span(x) 
with eigenvalue 0. But then Xa has distinct eigenvalues and so dim Oa = 3. □ 

Lemma 5.3. The maximum dimension of a proper face ofV is 6. 

Proof. Let F be a proper face of V and let S C SO (3) its stabilizer subgroup, 

S = {ge S0{3) I g.F C F} . 

This is a closed, proper subgroup, and thus either has dimension 1 (in which case it is a 
coaxial subgroup Qe or a rotation subgroup Q^), or it is finite and has dimension zero. 
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Let F° be the relative interior of F, those points of F which do not he in any other face 
of V of the same or smaUer dimension. If g.F° fl F° 7^ then g.F = F, and so g E S. 
Let dV be the boundary of V and consider the map /: 5*0(3) x F° — ^ dV defined by 

/ : {9,v) I — > g.v. 

This map is not 1-1: Suppose that g.v = h.w, for g,h E 5*0(3) and v,w E F°. Then 
h~^g.v = w and so h'~^g.F° fl -F° 7^ 0, which imphes that s := h^^g G 5*. Then s.v = w. 
This calculation shows that the fibers of / have the form 

{{gs-\ sv)\seS}, for g G S0{3) and x E F° . 

Thus we have the dimension calculation 

dimdV > dim^0(3) + dimF°-dim^. 

Since dimdV = 9 and dim 5*0(3) = 3, this gives 

(5.4) 6 + dim 5 > dimF. 

If 5* is finite, then dimF < 6. If 5* has dimension 1 so that it is either Qe or Q'^ for some 
e, then F could have dimension up to 7. By (12. 5p . W"^ = © © ?7| as a representation 
of 5*. If 5* = Qe, then Lemma [3171 implies that F has even dimension, and if 5* = Qjf, then 
Example 13.81 implies that F has dimension 0, 2, or 4, which completes the proof. □ 

5.1. Coaxial faces. A coaxial face of ^ is a face that is stabilized by some coaxial 
subgroup, Qe- By Lemma IX^ such a face is supported by a Qe-invariant linear function, 
which must factor through the projection to the trivial isotypical component of W"^, by 
Schur's Lemma. Since this component is (12. 5p . L is the puUback of a linear map 

3 (Ml, M2) I — > ttiMi + ttaMa E R . 

Up to a scalar, this is the composition of the Qe-invariant linear function Lg (14. ip on W 
with tTq,, which is the map L^^a defined by 

Suppose now that e is a unit vector. For each non-zero a G M^, define 
Mq := maximum eigenvalue of Xa , and 
rria := minimum eigenvalue of Xa ■ 
If X ^ V, then Xa ^ ^a, and so by Lemma 113] we have 

(5.5) M„ > Le,„(x) > 

with equality only when e is an eigenvector of Xa having eigenvalue Ma or ma- Thus 
coaxial faces are the faces of V defined by equality in (15.51) . 

For e G M'^ a unit vector and 7^ a G M^, define the coaxial face 

Fe,a := {X^V\ Le,,(x) = M„} . 

If X G Fe^a, then e is an eigenvector of Xa with eigenvector Mq. As in Section HJ each 
coaxial face F^^a is the convex hull of an orbit Qe-x'y for some x' ^ ^- 
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Theorem 5.6. Faces of V have dimension at most 6. The coaxial faces of V form a 
3 -dimensional family whose union is a 9-dimensional subset of the boundary of V if and 
only if Xi o-nd X2 have distinct eigenvectors. When this happens, almost all coaxial faces 
have dimension 6, have Caratheodory number 4, and are facets of V . 

If the boundary of V is the union of the coaxial faces, then Lemma 13.21 imphes that 
the Caratheodory number of V is at most 5, and we conjecture this is the case. If there 
are faces of dimension 6 that are not coaxial, then Caratheodory's Theorem implies that 
their Caratheodory number is at most 7. Then Lemma 13.21 implies the following corollary 
of Theorem 15.61 

Corollary 5.7. The Caratheodory number ofV is at most 8. 

By almost all in the statement of Theorem 15.61 we mean in the algebraic sense: Except 
for those a G lying in finitely many half-rays in M^, F^ a has dimension 6 when xi and 
X2 have distinct eigenvectors. The proof of Theorem 15.61 is done in the series of lemmas 
below. 

Remark 5.8. The condition that the magnetic susceptibility tensors Xi ^ind X2 have 
distinct eigenvectors has already been considered in protein folding. It implies that RDC 
measurements from the two ions are sufficient to remove the symmetry property of the 
RDC [131 . 

Since M_q, = —nia, there is no need for two types of coaxial faces as in Section HI Since 
if r > then M^a = rM^ and L^, = L_e, we have 

F = F = F 

^ e,a ^ —e,a ^ e,ra 

if r > 0. Thus we may assume that a lies on the unit circle in M^. We also only need 
to consider the unit vector e up to multiplication by ±1, that IS, cLS cL point in the real 
projective plane, MP^, which is a 2-dimensional manifold. 

Lemma 5.9. The coaxial faces Fe^a form a ?>- dimensional family parameterized by MP^ x 
. For each a & , any two coaxial faces F^^a O'nd Ff./ ^ are isomorphic. 

Since the boundary of V is 9-dimensional and it has a 3-dimensional family of coaxial 
faces, we see again that the maximum dimension of a coaxial face is 6. 
Proof. Suppose that e G is an eigenvector for Xa with maximal eigenvalue Mq,. Then 
X G Fe^a and F^^a is the convex hull of the orbit Qe-X- If -R G SO (3), then 

R.Qe.x = RQeR'.R.x = QRe.{R.x)- 

But {R.x)a is an anisotropic tensor having eigenvector Re with eigenvalue Mq,. Therefore 
FRe,a is the convex hull of Que-iR-x) = R-Qe-Xi and thus equals R.Fe^a- D 
We now determine the dimension of the coaxial faces. By Lemma 15. 9[ we need only 
study one coaxial face F^^a for each a E S^. We compute the dimension of the affine span 
of an orbit Qe-X^ where e is an eigenvector of Xa- This is the dimension of a coaxial face 
when the eigenvalue associated to e is a maximal eigenvalue of Xa- Since this dimension 
is the rank of a matrix, those entries are algebraic functions of a. Thus that for all but 
finitely many a, this rank will be constant and it will be smaller for a in that finite set. 
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Let a & and suppose that Xa,x' ^ Span(x) are linearly independent, and let e be 
a unit eigenvector of Xa- The decomposition fl2.5l) of W into Qe-iso typical components 
induces a decompostion of the tensors Xa and x' into their components in M © f/i © f/2, 

Xa = M„ © © t/„ and x' = M' ®x' ®y' . 

The f/i-component of Xa is 0, because e is an eigenvalue of Xa- Let di G {0, 1} be 
the dimension of the linear span of x' in Ui and 62 G {0, 1,2} be the dimension of the 
linear span of Uaill' in f^2- By Lemma 13.71 the dimension of the convex hull of Qe-X is 
di ■ dimf/i + d2 ■ dimf/2, which implies the following lemma. 

Lemma 5.10. The coaxial face Fe^a his dimension 2{di + c/2). 

Thus again a coaxial face has dimension at most 6. 

Lemma 5.11. Ifxi o-'ndx2 have a common eigenvector, then coaxial faces have dimension 
2 or 4. 

Proof. Fix a G S*^ and let / be a common eigenvector of xi and X2- Then it is an 
eigenvector of any Xa- Let x' be another tensor in Span(x) which is not proportional to 
Xa- Suppose that e = ei is an eigenvector of Xai that / G {61,62,63}, and write Xa and 
x' in the coordinates (12. 3p . 

Xa = 0,0, y,z) and x' = (M', w',x', y',z'). 

Note that {w',x') G Ui and {y,z), {y',z') G t/2. 

If / = 61, then {w', x') = (0, 0), and so di = 0. If / = 62 or 63, then z = z' = and so 
d2 = 1- In either case, c?i + ^2 < 3 and so the coaxial face Fg^^ has dimension 2 or 4. □ 

Lemma 5.12. If Xi (^''^d X2 do not have a common eigenvector, then there is a coaxial 
face with dimension 6. 

Proof. Suppose that 6 = 61 and 61, 62, 63 is an ordered basis of eigenvectors of Xi with the 
eigenvalue of e maximal. Write Xi and X2 in the coordinates (12.31) . 

Xi = x{Mu 0,0, 7,0) X2 = x{M2, w,x, y,z). 

The dimension of the coaxial face -Fei,(i,o) is 6 if and only if {w,x) 7^ (0,0) and 
(7, 0), (y, 2;) G U2 are linearly independent. Suppose that dimF < 6. We cannot have 
{w,x) = (0,0) for then 61 is a common eigenvector of Xi and X2, a contradiction. Thus 
the vectors (7, 0), {y, z) are dependent. 

If 7 = 0, then xi has a repeated smallest eigenvalue with eigenspace spanned by 62 and 
63. Changing the last two coordinates, we may assume that z = 0. We cannot also have 
either i;; = or x = for then Xi and X2 have either 62 or 63 as a common eigenvector. 
If ?/ = 0, then — X62 + we^ is a common eigenvector, so ?/ 7^ and -Fei,(i,o) has dimension 
4. In the coordinates (12.31) with respect to the ordered basis 63, 62, 61, — Xi and X2 are 

-Xi = X(^, 0,0, 3^,0) and ^2 = x(-^ - 2/, 0, a;, -3^ + |, w) . 

Since w,x 7^ 0, the affine span of Qe^-X has dimension 6. Since ^ is the maximal 
eigenvalue of — xi with eigenvector 63, this shows that Fe3,(_i,o) has dimension 6. 

The third possibility is that z = 0. But then the same arguments as in the previous 
paragraph show that coaxial face Fe3,(_i_o) has dimension 6. □ 



CONVEX HULLS OF ORBITS AND ORIENTATIONS OF A MOVING PROTEIN DOMAIN 19 



5.2. Structure and Caratheodory number of a coaxial facet. Suppose that F is a 
coaxial face of dimension 6. We may assume that F is the convex hull of the orbit Qe-X 
and that F spans the representation f/i © [/|. This Qe-orbit is the union of two orbits 
of its identity component Q^{— S0{2)). Call them and 0~ . By Example 13. 8[ each 
orbit spans a subrepresentation of Ui © f/| isomorphic to Ui © U2- Set 

F^ := convex hull of and := linear span of ~ Ui © f/2 • 

Proposition 5.13. The faces F^ each have dimension 4 and Caratheodory number 3. 
Points on their boundary are the convex hull of one or two vertices, while points in their 
relative interiors are the convex hull of three vertices. 

Proof. As in Example 13.81 we identify ~ 5*0(2) with the circle group 5^ and Ui,U2 
with C. Then z & acts on Ui as scalar multiplication by z and on U2 as scalar 
multiplication by z"^ and F^ has dimension 4. 

Let (m, f )^ G ~ t/i © ?72 be the point corresponding to x- Then 

0+ = {{e'%,e^''v) I < e < 27r} , 

and its convex hull is 

n n 

{ {Yl ^J-e'^' w, Yl ^J^^'^'^) I 5^ = 1' < ^1, . . . , ^„ < 27r} . 
But this is B.{u, v)"^, where B is the set of 2 x 2 diagonal matrices whose entries are 



(5.14) (J2 ' Y \^^'^') where = \ and < ^1, . . . , ^„ < 27 

i=i i=i 3 

Thus (and also F") is isomorphic to B. 

Curto and Fialkow [8] characterized the points of B. Let (a, 6) be a point of B (15.141) 
and p the corresponding measure on , 

»(e^^) = / ^ = ' 

1^ otherwise . 

Set 7jj := J^i z^z^dp{z) for < i,j with i + j < 2 and form the moment matrix 
M :: 



/700 


7oi 


710 ) 




( 1 


a 




710 


7ii 


720 


1 - 1 


a 


1 




\7oi 


702 


71 J 




[a 


b 





Proposition 5.15 (Curto and Fialkow [8]). The points {a,b) G B are exactly the points 
{a,b) G such that M is positive semi- definite. The rank of M is the minimum number 
of summands needed to represent the point {a,b) (15.141) . 

In particular, this implies that each body F^ has Caratheodory number 3. 
By Proposition 15. 15[ 

B = {{a,b)eC^\ l-|ap>0, 1 - |6p > 0, 1 + a% + a^b - 2\a\^ - \b\^ > O} . 

If 1 - |ap = then a e and so n = 1 in (I5.14p . If 1 - |6p = 0, then b e and either 
n = 1 in (15.141) or n = 2 with |6'i — 6'2| = vr. Thus if M has rank 3, then \a\ < 1 and |6| < 1, 
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and so (a, b) lies in the interior of B as the inequahties are strict. This imphes that points 
on the boundary of B are the convex hull of one or two vertices and this completes the 
proof of Proposition I5.13[ □ 

Now we complete the proof of Theorem 15.61 showing that the coaxial facet F has 
Caratheodory number 4. The coaxial face F is the convex hull of F~^ and F~ . Let v & F. 
We suppose that v ^ F^ U F~ , for otherwise v is the convex combination of at most three 
vertices. Then there exist G F^ and A G (0, 1) such that 

V = Xv^ + (1 - X)v' ■ 

If both f + and v~ lie on the boundary of their respective subfaces, then each is a convex 
combination of at most 2 vertices, and f is a convex combination of at most 4 vertices. 

Suppose instead that lies in the relative interior of F~^. The linear span of v and 
W~ has dimension 5 in the 6-dimensional space Ui © f/| and therefore it meets in 
a 3-dimensional affine subspace U^. Similarly the span of v and meets W~ in a 
3-dimensional subspace U~ . Observe that both and U~ contain the two-dimensional 
linear subspace fl W~ = Ui so that their span has dimension 4. 

Consider the cone over U~ fl F~ with vertex v. Removing v, this has two components. 
One meets F~ . Let C be the component which does not meet F~ , and let := Cfl . 
This is a convex set which contains and thus meets the relative interior of B~^ : = 
f/"*" n Points v' G fl B'^ are exactly those points of F'^ for which there exists a 
point v" G F~ such that v is a convex combination oiv' and v". There are two possibilities. 

(1) The boundary of C'^ meets the boundary of B~^. 

(2) Either the boundary of is a subset of B~^ or vice-versa. 

In the first case, let v' be a point common to the two boundaries. Then v' lies on 
the boundary of B^ and v" lies on the boundary of B~ . But these are subsets of the 
boundaries of F^, and so v is the convex combination of at most 4 vertices. 

In the second case, suppose that the boundary of B~^ is a subset of C~^. Since B~^ is 
the intersection F~^ with a hyperplane, its boundary must contain a vertex of F~^, as the 
set of vertices of F~^ is a connected 1-dimensional set whose convex hull is F~^. Suppose 
that v' G 5"*" is a vertex of F~^. Since v" G F~ is a convex combination of three vertices 
of F~ , we see that f is a convex combination of 1+3=4 vertices. If the boundary of 
is a subset of B^, then we may choose the point v' in the boundary of so that the 
corresponding point v" is a a vertex of F~. Again, v is the convex combination of 3 + 1 = 4 
vertices. □ 

We wish to thank Ivano Bertini, Claudio Luchinat and Giacomo Parigi of the Center 
for Magnetic Resonance of the University of Florence for suggesting and discussing with 
us this interesting problem. 
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